Multiplicative Noise Induces Zero Critical Frequency 
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Stochastic Bloch equations which model the fluorescence of two level molecules and atoms, NMR 
experiments and Josephson junctions are investigated to illustrate the profound effect of multiplica- 
tive noise on the critical frequency of a dynamical system. Using exact solutions and the cumulant 
expansion we find two main effects: (i) even very weak noise may double or triple the number of 
critical frequencies, which is related to an instability of the system and (ii) strong multiplicative 
noise may induce a non-trivial zero critical frequency thus wiping out the over-damped phase. 
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Many dissipative deterministic dynamical systems ex- 
hibit two phases of motion: an under-damped osciUa- 
tory behavior or an over-damped non-oscillatory motion. 
The transition between these common behaviors defines 
the critical frequency of the system, e.g. the critical fre- 
quency of the damped harmonic oscillator. Multiplica- 
tive noise is known to influence deterministic systems, 
in profound and surprising ways [l|, 13, 13, U, la, ISJ, |7J, ^ • 
Here we show how a stochastic perturbation induces a 
zero critical frequency for a particular non-trivial choice 
of noise strength, thus completely wiping out the over- 
damped phase. This might be counter intuitive at first 
glance, since we expect noise to work against oscillations, 
however as we soon demonstrate, in some cases the op- 
posite situation is found. The second interesting result 
we obtain is that even weak multiplicative noise may in- 
duce a doubling or a tripling of the number of critical 
frequencies of a system (in a way defined later) a result 
which is related to an instability of the noiseless dynam- 
ical system. Our results show how multiplicative noise 
may infiuence the critical frequency of a system in pro- 
found ways. 

We investigate the dynamics of the stochastic Bloch 
equation. The Bloch equation finds its applications, in 
many fields of Physics ranging from Nuclear Magnetic 
Resonance (NMR) [9|, llO| to single molecule spectroscopy 



11|, [lj| and Josephson's junctions [13|. We use the ex- 
ample of the optical Bloch equation, however with mi- 
nor modifications we may consider other systems e.g. 
magnetic systems. In particular we consider a two level 
electronic transition of an atom or a molecule, interact- 
ing with a continuous wave laser and a stochastic bath. 
The optical Bloch equation for Z[t) = [u, v, w, y), where 
(m, V, w) describes the usual Bloch vector, is 

|z(i) = M{t)Z{t) 



The initial condition is Z(0) = (0,0,-1/2,1/2) describ- 
ing a system in the ground state and y = 1/2 for all 
times. Here T is the radiate emission rate and 51 is the 
Rabi frequency describing the interaction of the tran- 
sition dipole of the system with the laser field (within 
rotating wave approximation). The stochastic detuning 
5L{t) — ujl — (wo + v{t)) describes the interaction of the 
system with the bath in the spirit of the Kubo- Anderson 
line shape theory [ll|, [ij, ll5| , namely 5l (t) is a stochas- 
tic process describing spectral diffusion, uj^ is the laser 
frequency and loq is the absorption frequency of the two 
level system. Spectral diffusion is found in many molecu- 
lar, atomic and magnetic systems and is well investigated 
[111 . Ilq . The noise is called multiplicative since in Eq. 
(fT|) the spectral diffusion process multiplies the vector 
Z. Beyond spectral diffusion the equations describe a 
two level system in the process of resonance ffuorescence, 
where the laser frequency exhibits fluctuations, namely 
the detuning is a random function of time. 

For the noiseless case ^{t) — and for zero laser de- 
tuning Wi = ujq we have a simple damped harmonic os- 
cillator for w 
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When ri is larger than the critical frequency ilc = r/4 
the system exhibits under-damped Rabi oscillations while 
when rj < J7c it decays to the steady state monotonically. 
The critical frequency provides the quickest approach of 
the amplitude of the damped harmonic oscillator to zero. 

Now we consider the system in the presence of the 
spectral noise and investigate the average behavior (w). 
What will happen to the critical frequency of the system? 
and can we choose parameters of the noise in such a way 
that the critical frequency of the noisy system is zero? 

The formal solution of the problem is given in terms 
of the time ordered exponential 



M(i) = 



f -^ 5L{t) 

-&L{t) -^ -Vi 

VI -V -T 

V 



(1) 



im) 



T 



exp 



M(T)dT]\)m. (3) 



In practice it is generally difficult to find explicit solutions 
due to the combination of the time ordering operator T 
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FIG. 1: The phase diagram for the critical frequencies of the 
optical Bloch equation with multiplicative two state noise. 
The 0,1,2,3 indicate the number of critical frequencies, as 
defined in the text. In (a) we see a line of zero critical fre- 
quency in the fast modulation limit v < R, which is well 
approximated by the cumulant expansion (the dashed line 
v/V = yj R/V). In the absence of spectral diffusion z/ = we 
have a single critical point, thus as shown in (b) the addition 
of weak noise v/V << 1 may strongly influence the critical 
frequency in the sense that we find there a phase with 2 or 3 
critical frequencies. 

and the average over the multiplicative stochastic pro- 
cess denoted w\\h (• • •) in Eq. ([3]). Here we find an exact 
solution for a dichotomic two state Kubo- Anderson pro- 
cess ilSJ]. With this solution we will explore whether the 
motion is over-damped or under-damped. We later show 
that our findings are general beyond the exactly solvable 
two state process. In particular we consider i'{t) = i'h{t), 
where h{t) = +1 or h{t) = — 1 describes the stochastic 
two state process with a rate R for transitions between 
-|-1 and —1. Such a model is applicable in single molecule 
spectroscopy in glasses [12|, |l7|, |l8[ and was used exten- 
sively for the line shape theory of Kubo and Anderson. 

We use Burshtein's method [19] of marginal averages 
to solve the Kubo-Anderson process with zero laser de- 
tuning lol = wq. We define the marginal average vector 
{{Z{t))+, {Z{t))-) which is the average of Z{t) given that 
at time t the stochastic process had the value ±1 corre- 
spondingly. The equation of motion for the marginal 
average vector is 



d_ ( {Z) + 
dt \ (Z>_ 
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where / is the identity matrix. The operators Aq and 
Aq are Bloch matrices corresponding to the state of the 
spectral diffusion Sl = +1^ or Sl — —v respectively. To 
solve the problem we must diagonalize the 8x8 matrix in 
Eq. (j4]), then complex eigenvalues yield under-damped 
oscillatory modes while real eigenvalues correspond to 
over-damped modes. The eigenvalues {A} are found us- 
ing the characteristic polynomial of Eq. ([4]) 



A(A + 2i?)Pi(A)P2(A)=0, 



where two cubic polynomials are defined as 



(6) 



Pi(A)= (8r!2 + 4(r + A)(r + 
+2(r 
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P2(A)= 8(r + 2A)i?2+ 

+2 (4z^2 + (r + 2A)(3r + 4A)) R+ 

+2(r + 2A)r!2 + (r + A) ((r + 2A)2 + Av'^) . 

(7) 
We have thus reduced the problem to finding the roots 
of two third order polynomials -Pi (A) = and -P2(A) = 0. 
After some algebra one can show that only the roots of 
Pi (A) enter the solution of {Z) = (Z)^ + {Z)_ as well 
as the eigenvalue A = [see Eq. ([5])] which yields the 
steady state solution [2C1J . 

A physical observable is the intensity of emitted light 
(I(t)) which is equal to F times the population in the 
excited state {I (t)) = T{^{w) + ^). The eigenvalues 
{Ai,A2,A3} are the solutions of Pi (A) = and the in- 
tensity is 
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(8) 
with FsD = v^ JR. The steady state solution is 
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This expression when fi — > is the well known Kubo- 
Anderson line shape at zero laser detuning. 

We now focus our attention on the eigenvalues {A.^} to 
determine whether the solution is over-damped or under- 
damped. The motion is called over-damped if all eigen- 
values {\i\ are real otherwise it is under-damped. The 
condition for over-damped behavior is that the discrimi- 
nant T> of Pi (A) be less than zero, explicitly we have 
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Recall that for the noiseless case we have a single crit- 
ical frequency flc — r/4. As shown in Fig. [T] in the 
presence of multiplicative noise the phase diagram of the 
motion is very rich: 

(i) In the slow modulation v ^ R strong coupling 
V /T ^ 1 regime, the solutions are always oscillatory and 
the critical frequency is (indicated by in Fig[T^). 
(ii) When v Ki R we obtain over-damped motion when 
fici < ^ < ^C2 so we have two critical frequencies (in- 
dicated by 2 in Fig. [1^). 

(iii) In the fast modulation R > v strong coupling 
v/T > 1 limit, we find a single critical frequency simi- 
lar to the noiseless case, except for a surprising line on 
which Q.C = (denoted with in Fig. [T^). 
(iv) In the weak noise limit v <^T, the solution yields 
either 2 or 3 critical frequencies (see Fig. [TJd) . Thus weak 
noise modifies the solution dramatically by doubling or 
tripling the number of critical frequencies of the system. 

To understand better the phase diagram Fig. [T] we 
present in Fig. [2] the behavior of the solution of the op- 
tical Bloch equation in the absence of the multiplicative 
noise i.e. u = Q. The figure is a phase diagram in the Rabi 
frequency Jl/F and detuning {uol — ^)/^ plane showing 
the regions of over-damped and under-damped behavior. 
Fig. [5] illustrates that the noiseless system is unstable 
in the sense that for any small detuning and low enough 
Rabi frequency we get an oscillatory behavior while for 
zero detuning the solution is over-damped. This insta- 
bility of the noiseless solution explains why even adding 
a weak perturbation v -^T strongly affects the system. 
Namely, for weak noise (Fig, [TJd) wc find either 2 or 3 
critical frequencies instead of 1 for the noiseless case. 

As mentioned before for slow modulation v > R and 
strong coupling ^' ^ F the motion is always under- 
damped. To understand this behavior we again refer 
to the noiseless case presented in Fig. [2l where we ob- 
serve that large detuning means an oscillatory solution. 
Namely oscillations for the noise free Bloch equation are 
induced by two mechanisms, the Rabi frequency and the 
detuning. Hence it is not surprising that strong and slow 
noise in Fig. [Iji (i.e. i/ > R,T) may induce oscillations 
and the wipe out of the over-damped motion. 

Far less trivial is the wipe out of over-damped motion 
in the fast modulation limit, i.e. the line of zero critical 
frequency in Fig. [IJi. To investigate this behavior we 
consider the limit i? — > oo. Then using Eq. (jlOp we find 
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FIG. 2: Phase diagram of the optical Bloch equation in the 
absence of the multiplicative noise. The darker area is the 
over-damped phase. For zero detuning ujl~^o = the critical 
frequency is fie ~ r/4. Notice the cusp at zero detuning 
which makes the solutions unstable to multiplicative noise. 



the critical frequency 
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where the limit is taken with FgD remaining finite. We 
see that £7^ = when Fgo = F, namely when v/T = 
y^R/r. This line is shown in Fig. [l]as a dashed line. 

Expanding the exact solution in Fsd, one can show 
that for R > F/8 any amount of noise will lead to a 
decrease of ftc according to 
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where the leading F/4 term describes the noiseless case. 
The decrease of Jlc is explained by the fact that the noise 
removes the system from zero detuning and hence solu- 
tions tend to be more oscillatory (i.e. the critical fre- 
quency is reduced). The surprising result is that by in- 
creasing the noise level we reach a limit where the critical 
frequency is zero. Such a behavior in the fast modulation 
limit could not be anticipated without our mathematical 
analysis. The behavior of fie is illustrated in Fig. [3l 
which shows the decrease of the critical frequency until 
it reaches the value J7c = 0. 

It is natural to ask if the behavior we found is general 
or limited to the example of a two state process. For this 
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FIG. 3: The critical frequency flc as a function of Fsd which 
is the measure of noise strength. For small noise levels, i.e. 
Fsd/F < 1, the critical frequency flc decreases as anticipated 
in Eq. (|12|) . The figure illustrates the existence oi Qc ~ 
for a particular noise value. We show the critical frequency 
obtained from the exact solution (dot-dashed) and the cumu- 
lant approximation (dashed) for R — 2.5r. The solid line is 
the critical frequency f2c at _R -^ oo Eq. (Ill|l . 



aim we have used the cumulant expansion 21|, to inves- 
tigate the critical frequency flc of the system. We con- 
sider a stationary process h(t) whose correlation function 
is {h{t)h{t + t)) = exp(— i?r). The cumulant expansion 
works well when the Kubo number f/R is small. Within 
this approximation 21^ 
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where Aq — A^ \ ^^q and 
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with Ai = 



Solving the 



cr \ , / 1 

oj ^^^^= (,-1 
integrals leads to cumbersome equations for 51c- However 
in the limit i? — > oo and z^ ^ cxd in such a way that 
TsD = v^ /R remains constant we find that Eq. pT|) is 
valid and therefore the equation is not limited to the two 
state Kubo- Andersen model. This means that for a large 
class of stochastic processes, multiplicative noise induces 
zero critical frequency in the fast modulation limit and 
the curve v /T = y^R/T, on which flc — shown in Fig. 
[ija), is a general behavior. 

To further validate the generality of our results we have 
solved semi-analytically and with the help of Mathemat- 
ica: (i) two state model with two non identical rates de- 
scribing the transitions between up and down states, and 
(ii) models with three states. These models show behav- 
iors similar to our findings. 

The same effects cannot be found for linear systems 
driven by additive noise with zero mean, since the aver- 
aged equations have the same critical frequency as the 



noiseless case. Hence for linear systems the effects we 
have found are limited to systems with multiplicative 
noise. However nonlinear systems with additive noise 
may exhibit behaviors identical to those investigated in 
this manuscript. To see this we add a coordinate 5 to the 
description of the system thus Eq. ([T|) is written as: 



d_ 
dt 
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This equation is a non linear stochastic equation with 
additive noise (i.e. ^(i) in Eq. (fTSJ) ) which is equiva- 
lent to the multiplicative Eq. ([T]). We see that the two 
main effects found in this manuscript: (a) noise induc- 
ing zero critical frequency and (b) the doubling or the 
tripling of the number of critical frequencies, even for 
weak noise, may be found either in linear multiplicative 
systems or non linear additive systems. Thus we expect 
the main features of our results to be valid for a vast class 
of stochastic dynamical systems. 
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